Sequential generation of entangled multi-qubit states 
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We consider the deterministic generation of entangled multi-qubit states by the sequential coupling 
of an ancillary system to initially uncorrelated qubits. We characterize all achievable states in terms 
of classes of matrix product states and give a recipe for the generation on demand of any multi-qubit 
state. The proposed methods are suitable for any sequential generation-scheme, though we focus on 
streams of single photon time-bin qubits emitted by an atom coupled to an optical cavity. We show, 
in particular, how to generate familiar quantum information states such as W, GHZ, and cluster 
states, within such a framework. 
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Entangled multi-qubit states are a valuable resource 
for the implementation of quantum computation and 
quantum communication protocols, like distributed 
quantum computing Q, quantum cryptography 2] or 
quantum secret sharing Using photonic degrees of 
freedom as qubits, say, polarization states or time-bins 
of energy eigenstates, has the advantage that photons 
propagate safely over long distances. Consequently, pho- 
tonic devices are the most promising systems for quan- 
tum communication tasks. For this purpose, a lot of 
effort has been made in recent years to develop efficient 
and deterministic single photon sources 0, H, S IS • 

Photonic multi-qubit states can be generated by let- 
ting a source emit photonic qubits in a sequential manner 
|lO] . If we do not initialize the source after each step, the 
created qubits will be in general entangled. Moreover, if 
we allow for specific operations inside the source before 
each photon emission, we will be able to create different 
multi-qubit states at the output. In fact, this is a partic- 
ular instance of a general sequential generation scheme, 
where an ancillary system is coupled in turn to a number 
of initially uncorrelated qubits. 

It is the purpose of this paper to provide a com- 
plete characterization of all multipartite quantum states 
achievable within a sequential generation scheme. It 
turns out that the classes of states attainable with in- 
creasing resources are exactly given by the hierarchy of 
so-called matrix-product states (MPS) 0, E^. These 
states typicall y a ppear in the theory of one-dimensional 
spin systems 'l3], as they are the variational set over 
which Density Matrix Renormalization Group techniques 
are carried out Thus, our analysis stresses the im- 
portance of MPS, since we show that they naturally ap- 
pear in a completely different and relevant physical con- 
text. Moreover, particular instances of low-dimensional 
MPS, like cluster states 15] or GHZ states 1 161 . are a 
valuable resource in quantum information [13. Con- 
versely, we will provide a recipe for the generation on 
demand of any multi-qubit state within a sequential gen- 



eration scheme. Due to the general validity of these re- 
sults, we will first state and prove them without refer- 
ring to any particular physical system. This will be then 
applicable to all sequential setups, like streams of pho- 
tonic qubits emitted either by a cavity QED (CQED) 
source 0, 0, ^ or by a quantum dot coupled to a mi- 
crocavity [a^lfll- 

In the second part, we will focus on the physical im- 
plementation of these ideas within the realm of CQED. 
The role of the ancillary system will be performed by a 
D-\eYe\ atom coupled to a single mode of an optical cav- 
ity. The sequentially generated qubits will be time-bin 
qubits |0) and |1), describing the absence and presence 
of a photon emitted from the cavity in a certain time 
interval (see Fig. 1). 
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FIG. 1: A trapped D-level atom is coupled to a cavity qubit, 
determined by the energy eigenstates |0) and |1). After ar- 
bitrary bipartite source-qubit operations, photonic time-bins 
are sequentially and coherently emitted at the cavity output, 
creating a desired entangled multi-qubit stream. 

We will concentrate on setups where all intermediate 
operations are arbitrary unitaries and the ancilla decou- 
ples in the last step. The latter enables us to generate 
pure entangled states in a deterministic manner without 
the need of measurements. Let Ha — and He — 
be the Hilbert spaces characterizing a D-dimensional an- 
cillary system and a single qubit (e.g. a time-bin qubit) 
respectively. In every step of the sequential generation of 
a multi-qubit state, we consider a unitary time evolution 
of the joint system Ha ^ ^B- Assuming that each qubit 
is initially in the state |0) (i.e., the time-bin is empty), we 
disregard the qubit at the input and write the evolution in 
the form of an isometry V : Ha Ha^Hb- Expressing 
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the latter in terms of a basis V = ^ ^ ^a,/3l^' 
isometry condition reads X]i=o = 1, where each V* 

is a I) X I) matrix. We choose a basis where 
are any of the D ancihary levels. If we now apply suc- 
cessively n, not necessarily identical, operations of this 
form to an initial state \^i) ^ we obtain the state 

l*)=l^[n]...l^[2]V[l]|^/). (1) 

Here, and in the following, indices in squared brackets 
represent the steps in the generation sequence. The n 
generated qubits are in general entangled with the ancilla 
as well as among themselves. Assuming that in the last 
step the ancilla decouples from the system, such that 
1^) = I^f) ^ IV^), we are left with the n-qubit state 

1 

l^>= E {VF\V{-y..V{^^\^l)\in,...,il). (2) 

ii...in=0 

States of this form are called matrix-product states 

(MPS) HI 113, 

and play a crucial role in the theory of 
one-dimensional spin systems. Equation ^ shows that 
all sequentially generated multi-qubit states, arising from 
a D-dimensional ancillary system Ha^ are instances of 
MPS with D X D matrices and open boundary condi- 
tions specified by \(pi) and {(pr)- We will now prove that 
the converse is also true, i.e., that every MPS of the form 

\rl,) = {^F\V[n]...V[,0i), (3) 

with arbitrary maps V[fe] : Ha Ha ^ Hb^ can be gen- 
erated by isometrics of the same dimension, and such 
that the ancillary system decouples in the last step. 
Since every state has a MPS representation, this is at 
the same time a general recipe for its sequential genera- 
tion. The idea of the proof is an explicit construction of 
all involved isometrics by subsequent application of sin- 
gular value decompositions (SVD). We start by writing 
{^F\V[n] = V^^|M[^], where the 2x2 matrix V^^^ is the left 
unitary in the SVD and M[^] is the remaining part. The 
recipe for constructing the isometrics is the induction 

(M[fc] ® h)V[k-i] = V['fc_i]M[fe_i], (4) 

where the isometry ^[/^-i] is constructed from the SVD of 
the left hand side, and M[k-i] is always chosen to be the 
remaining part. After n applications of Eq. in Eq. 
from left to right, we set 1^:'/) = M[i]|(^/), producing 

m = v{„]...v{,^\ipi). (5) 

Simple rank considerations show that has dimen- 

sion 2min[D,2^] x min[D,2^+^]. In particular, every 
V^'^l could be embedded into an isometry V^^] of dimen- 
sion 2D X D. Physically, this just means we would have 



redundant ancillary levels that we need not to use. Fi- 
nally, decoupling the ancilla in the last step is guaranteed 
by the fact that, after the application of Vj^_i], merely 
two levels of Ha are yet occupied, and can be mapped 
entirely onto the system Hb- This is precisely the action 
of V[n] through its embedded unitary V^^^ . 

This proves the equivalence of two sets of n-qubit 
states, which are described either as D-dimensional MPS 
with open boundary conditions, or as states that are gen- 
erated sequentially and isometrically via a D-dimensional 
ancillary system which decouples in the last step. Mo- 
tivated by current cavity QED setups, we will now pro- 
vide a third equivalent characterization, namely, a set of 
multi-qubit states that are sequentially generated by a 
source consisting of a 2D-level atom. In contrast to the 
first sequential scheme, the latter will not require arbi- 
trary isometrics. 

Consider an atomic system with D states \ai) and D 
states \bi), so that Ha = Ha^Hb ^ C^(8)C^ That is, we 
will write \(p)\l) for a superposition of |ai) states, whereas 
\(p)\0) corresponds to a superposition of \bi) states. Since 
the last qubit marks the atomic level, whether it belongs 
to the \ai) or to the \bi) subspace, we will refer to it as the 
tag-qubit and write Ha = Ha' ^ Hr- Now consider the 
atomic transitions from each \ai) state to its respective 
\bi) state accompanied by the generation of a photon in 
a certain time-bin. This is described by a unitary evolu- 
tion, since now called "D-standard map" , of the form 

T: \^)a'\1)t\0)b ^ \^)a'\0)t\1)b , 

\^)a'\0)t\0)b ^ \^)a'\0)t\0)b ■ (6) 

Hence, T effectively interchanges the tag-qubit with the 
time-bin qubit. If, additionally, arbitrary atomic uni- 
taries Ua are allowed at any time, we can exploit the 
swap caused by T in order to generate the operation 

V\^) = (0|rT(c/^(|(p)^,|0)r)|0)H) , (7) 

which is the most general isometry V : Ha' Ha' ^ 
Hb' Therefore, the so generated n-qubit states include 
all possible states arising from subsequent applications of 
2D X D-dimensional isometrics. On the other hand, they 
are a subset of the MPS in Eq. (O with arbitrary 2D x D- 
dimensional maps, assuming that the atom decouples at 
the end. Hence, these three sets are all equivalent. 

Now, we show how these results can be applied in the 
realm of cavity QED, where an atom is trapped inside a 
high-Q optical cavity, and we aim at generating multi- 
photon entangled states. A laser may excite the atom, 
producing subsequently a photon in the cavity mode, 
which, after some time, is emitted outside the cavity (Fig. 
1). We consider two different scenarios, corresponding to 
the two families of states considered above. First, we may 
have fast and complete access to the atom-cavity system. 
In consequence, after the implementation of the desired 



isometry in each step, we should wait until the photon 
leaks out of the cavity before starting the next step. In 
this case, according to the analysis above, we will be able 
to produce arbitrary D-dimensional MPS with D equal 
to the number of involved atomic levels. Second, we may 
have a 2D-level atom {D \ai) levels and D \hi) levels) 
and two kind of operations: (i) fast arbitrary operations 
which allow us to manipulate at will all atomic levels; 
(ii) an operation which maps each \ai) state to its corre- 
sponding \hi) state while creating a single cavity photon, 
allowing a taylored output. Here, we will also be able to 
produce arbitrary D-dimensional MPS. 

In the following, we will illustrate the above statements 
with a specific example which is based on present cav- 
ity QED experiments [1, 0, 0- We consider a three- 
level atom coupled to a single cavity mode in the strong- 
coupling regime. An external laser field drives the tran- 
sition from level \a) to the upper level \u) with coupling 
strength (^o, and the cavity mode drives the transition 
between \u) and level \h) with coupling strength ^f, in a 
typical A configuration [see Fig. 2(a)]. We choose the 
detunings A, with |A| ^ {g^Vto}^ and assume that the 
cavity decay rate n is smaller than any other frequency in 
the problem, so that we can ignore cavity damping dur- 
ing the atom-cavity manipulations. By eliminating level 
\u)^ we remain with an effective D = 2 atomic system 
plus the cavity mode. We will show how, by control- 
ling the laser frequency and intensity, it is possible to 
generate arbitrary 2-dimensional MPS. Note that, by al- 
lowing the manipulation of D effective atomic levels, it is 
straightforward to extend these results to the generation 
of D-dimensional MPS. 

According to the results presented above, we just have 
to show that we can implement any isometry V : Ha 
Ha^Hb- In fact, we will show how it is possible to imple- 
ment arbitrary operations on the atomic qubit and the 
V SWAP operation on the atom-cavity system, which 
suffice to generate any isometry V (since they give rise 
to a universal set of gates for the two qubit system 18]). 
The atomic qubit can be manipulated at will using a 
Raman laser system, as it is normally done with trapped 



ions [l^ |2^ . In order to implement the ^/ SWAP, we no- 
tice that the atom-cavity coupling is described in terms of 
the Jaynes-Cummings model [see Fig. 2(b)], where the 
coupling constant is controlled by the laser. Thus, 
application of laser pulses with the appropriate duration 
and phase [23, will implement the unitary operation 
U = e"'^, where generator G = ( |a,Q)(6, 1| + H.c.)7r/4, 
which corresponds to the desired V SWAP operation. In 
order to generalize this scheme to an arbitrary D-level 
system, we notice that we can view the atom as a set 
of M qubits (with D < 2^). Thus, if we are able to 
perform arbitrary atomic operations, together with the 
V SWAP operation on two specific atomic levels as ex- 
plained above, we can then implement a universal set of 
gates and, in consequence, any arbitrary isometry. 




FIG. 2: (a) Atomic level structure: levels \a) {\b)) and \u) 
are coupled by a laser (cavity mode) off resonance, (b) After 
adiabatic elimination of the upper state \u), we are left with 
a Jaynes-Cummings type of Hamiltonian, where states |a,n) 
and \b,n + 1) are coupled. Both, the energy difference of 
those levels and the corresponding Rabi frequency depends 
on n. The reason for the first is the AC-Stark shift, whereas 
the second is due to the Jaynes-Cummings coupling. 



In the rest of the paper, we will use another setup 
which is closely related to current experiments 0, 0, H 
and optimizes our second method for MPS generation. 
In this frame, we will show how to gene rate familiar 
multi-qubit states like W 22], GHZ [wj . a nd cluster 
states 15], which are all MPS with L> = 2 0. 

For the purpose, we consider an atom with three effec- 
tive levels I62)} trapped inside an optical cavity. 
With the help of a laser beam, state \a) is mapped to the 
internal state and a photon is generated, whereas the 
other states remain unchanged. This physical process is 
described by the map 



(8) 



Mab: \a) ^ \bi)\l) 
\bi) ^ \bi)\0) 
I62) ^ |62)|0) 



and can be realized with the techniques used in 
After the application of this process, an arbitrary opera- 
tion is applied to the atom, which can be performed by 
using Raman lasers. The photonic states that are gener- 
ated after several applications are those MPS where the 
isometrics are given by V^^ = MabU^ , with i = 1, . . . , n, 

Uj^ being arbitrary unitary atomic operators. 

For example, to generate a W-type state of the form 

\^w) = e^*^^sin6i|0...01) + coseie^*^2sin62|0...010) 

+ cos6i...cos6n-2e'*^"-' sin6n-i|010...0) 
+ cos6i...cos6n-i|10...0), (9) 

we choose the initial atomic state \(pi) = I62) and opera- 
tions U]^ = U^l^i^i, Oi), with z = 1, . . . , n - 1, where 

U]S{^i.&i) = cosei\k){k\ ^ cosei\l){l\ ^ e'^^smei\k){l\ 



— e 



sin6i|/)(A:| + \m){m\ 



(10) 



and {/c,/,m} = {a, 61, 62}- To decouple the atom from 
the photon state, we choose the last atomic operation 
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^a'^ = ^ah2^^^^/^) a^d' af^^^ the last map M^^, the 
decoupled atom will be in state \hi). 

To produce a GHZ-type state in similar way, we choose 
1^,) = |a), f/W = C/X(*i,ei), U^^ = uX{0,^m, with 
i = 2, . . . , n - 1, and = ,^(0, 7r/2)[/,^g^(0, 7r/2). 

For generating cluster states, we choose \(pi) = \b2), 
U^A = UXi^i, OiWX (0' ^/2)' with i = 1, . . . , n-1, and 
U^A = Kl ^n)U§^b, (0, 7r/2)UX (0, 7r/2), obtaining 

n 

IV') = (g)(o°-i|o), + ol_i|i)0, (11) 

where 0?_^ = cos 6i|0)i_i(0| - e"^*^^ sin 6i|l)i_i(l| 
and 0^_i = e'^^s'mQi\0)i-i{0\ + cosei|l)i_i(l|, with 
i = 2,...,n — 1. Operators 0^_i and O^Li act on the 
nearest neighbor-qubit i — 1 under the assumption Oq = 
cos©i and Oq = ^^^^ sin 61. If one chooses = and 
6i = 7r/4 this leads to the cluster states defined by 

1 

IV^ci) = ^ K-ilO). + , With ^ 1. (12) 

The formalism presented here is also valid for other 
types of single-photon sources, in the context of cavity 
QED or quantum dots. For example, it could be ex- 
tended to characterize the polarization-entangled multi- 
qubit photon states generated by an analogous cavity 
QED photon source 10]. In fact, the presented ideas 
and proofs apply to any multi-qudit state with Hb — 
that is generated sequentially by a /^-dimensional source. 

In a wider scope, we have established a formalism de- 
scribing a general sequential quantum factory, where the 
source is able to perform arbitrary unitary source-qudit 
operations before each qudit leaves. Apart from the mul- 
tiphoton states, the present formalism applies also to 
many other physical scenarios: (a) to coherent microwave 
cavity QED experiments 23], where atoms sequentially 
cross a cavity, and thus the outcoming atoms end up in 
a MPS with the dimensions given by the effective num- 
ber of states used in the cavity mode; (b) a light pulse 
crossing several atomic ensembles 12 4| . which will be left 
in a matrix product Gaussian state 25]; (c) trapped ion 
experiments where each ion interacts sequentially with 
a collective mode of the motion 19, 2^|2^. Note also 
that one can include dissipation in the present formalism, 
by replacing; MPS by matrix product density operators 
I12L l27j . This description applies, for example, to the 
micromaser setup 28] and other realistic scenarios. 
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